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1. Introduction model atmosphere. The unceitainty in parall ax and bolometric 

in flux do minate the errors. As a recent example. iNeilson & Lesteil 

O The basic properties and evolutionary status of a star can be ^^m) determined the fundamental parameters of three M- 

^ determined knowing its mass, luminosity, radius and chemical ^yp^ ^^^^ j^^^ combined infrared interferometry with spectro- 

T-H , composition. Nowadays we are reaching a pomt where mass photometric observations. The interferometric data were ob- 

^ . estimates are available not only for binary but alsofor single ^^^^^^ ^j^^ the VLTl/VlNCl instmment, and previously in- 

k> objects (see the recent review bylAufdenbergetal., 2009). This terpreted by using Phoenix an d Atlas atmospheric models 

major advance has been mainly achieved by establishingjnter- ,'Hauschildt e t alL^ 'KuruczV 1993: respectively) in the work 
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ABSTRACT 

Context. Giant stars, and especially C-rich giants, contribute significantly to the chemical enrichment of galaxies. The determination 
of precise parameters for these stars is a necessary prerequisite for a proper implementation of this evolutionary phase in the models of 
galaxies. Infrared interferometry opened new horizons in the study of the stellar parameters of giant stars, and provided new important 
constraints for the atmospheric and evolutionary models. 

Aims. We aim to determine which stellar parameters can be constrained by using infrared interferometry and spectroscopy, in the case 
of C-stars what is the precision which can be achieved and what are the limitations. 

Methods. For this purpose we obtained new infrared spectra and combined them with unpublished interferometric measurements for 
five mildly variable carbon-rich asymptotic giant branch stars. The observations were compared with a large grid of hydrostatic model 
atmospheres and with new isochrones which include the predictions of the thermally pulsing phase. 

Results. For the very first time we are able to reproduce spectra in the range between 0.9 and 4 yum, and K broad band interferometry 
with hydrostatic model atmospheres. Temperature, mass, log(g), C/0 and a reasonable range for the distance were derived for all 
the objects of our study. All our targets have at least one combination of best-fitting parameters which lays in the region of the HR- 
diagram where C-stars are predicted. 

Conclusions. We confirm that low resolution spectroscopy is not sensitive to the mass and log(g) determination. For hydrostatic 
objects the 3yum feature is very sensitive to temperature variations therefore it is a very powerful tool for accurate temperature 
determinations. Interferometry can constrain mass, radius and log(g) but a distance has to be assumed. The large uncertainty in the 
distance measurements available for C-rich stars remains a major problem. 

Key words, stars: AGB and post-AGB - stars: atmospheres - stars: carbon - stars: fundamental parameters - techniques: spectroscopic 
- techniques: interferometric 



H • ferometry as a standard tool for investigating stars (Wittkowskil Wittkowskiltdr(l2004, 2006a b). Neilson & Lester (20o|) 

.y: 12004). The combination of high angular and spectral resolution ^^^^ ^^e new generation of S Atlas models (ILester & Neilsc^ 

gives the possibility to study the spatial structure of single stellai- 20O8) for the interpretation of the data: in this way the authors 

objects, which were before ti-eated only as point sources. Red gi- ^^^^ ^^j^ ^^jy determine the fundamental parameters, but 

ant stars are very good targets for interferometric investigations ^j^^ constrain the atmospheric models used. The typical ap- 

because of their extended atmosphere and their brightness in the p^.^^^^ followed in these works to determine the stellar pai-am- 

infrared. These stars are main contnbutors to the infrared light ^^^^^ j^. determination of the limb darkened radius by fitting 

and to the chemical enrichment of galaxies. Therefore an accu- ^^e interferometric observations, the resulting radius is converted 

rate determination of their fundamental parameters is mandatory Rosseland radius; (ii) a linear radius is derived by assuming 

for a proper implementation in the models of galaxies. ^ distance; (iii) an efi^ective temperature is derived by fitting the 

I Wittkowski et al.| (12001|) made a first attempt to derive the pa- ^ED; (iv) the luminosity is obtained using radius and T.ff, (v) the 

rameters of a sample of asymptotic giant branch (AGB) stars, by gj-^^jty ^^d the mass of the star are constrained by comparison 

comparing the spectro-interferometric observations with Kurucz ^j^j^ jj^e theoretical models 
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The targets of our study are an important subclass of the 
AGB objects: the carbon-rich AGB stars. According to theory 
only stars with a range of mass between 1 and 4 Mq undergo 
the third dredge-up during their AGB phase, with the result of 
increasing their C/O ratio ( Iben & Renzinil Il983l) . The spectra 
of these giants are consequently dominated by the absorption 
features of carbonaceous molecular species like: CO, CN , HCN, 
C2, C 3, C2H2 (Tsuii, 1986; J0rgensen et al., 2000; Lo idl et al.L 
I2OOI . As a result, these objects are main contributors to the C- 
enrichment of the ISM. 

The state-of-the-art model atmospheres for C-stars are 
very pro mising, being abl e to produce la rge grids of models 
dAringer et al., 2009; Matts son etal.L 12010 ) which describe the 
observed properties of these objects reasonably well. The hydro- 
static model atmospheres for C-stars were c ompared with optical 
and infrared spectroscopy in lj0r^ensen et al. (2000); Loi dl et aP 
(l200ll) ; lAbia et al.1 (12010 ) and references therein. The major dif- 
ference between an earlier generation o f models as used by 
(|j0rgensen et al.Ll2000HLoidl et al.Ll200lh and the hy drostatic at- 
mospheres applied in our work ("Arin ger et al.l l2009l ) is the inclu- 
sion of atomic opacities and an update to new and more accurate 
molecular opacities. 

The stellar evolutionary calculations reached a highly so- 
phisticated level as well, being able to include a detailed ther- 
mally pulsing phase with third dredge-up, hot-bottom burning 
and variable molecular opacity (Marigo et al., 2008). 

The purpose of this work is to find out which stellar param- 
eters can be determined for weakly variable C-stars (visual am- 
plitude < 2 mag) by using spectroscopy, infrared interferometry 
and hydrostatic model atmospheres; and what is the accuracy. 
We aim to test the stellar atmospheric models with different tech- 
niques and a multiwavelength approach. This is a very important 
step since errors in the model structure limit quantitative checks 
of stellar evolution. Our approach in the paramet er determination 
is sl ightly different from the one adopted by Neilson & Lester! 
(12008) for a few simple reasons: (1) the Rosseland radius is not 
a direct observable, and in the C-stars (to the contrary of M- 
type AGB stars) there are no windows in t he near i nfrared where 
to measure a "continuum" radius (Paladi ni et al.l |2009); (2) it 
is very rare to obtain simultaneous multiwavelength photomet- 
ric observations, which is mandatory when dealing with variable 
stars. Nevertheless the parameters determined in this work will 
be compared in a follow-up paper (van Belle et al., in prep.) with 
the one obtained by combining photometric and interferometric 
observations. 

A short description of the observations acquired and of the 
data reduction for both spectroscopy and interferometry will 
be given in Sect. |2l The hydrostatic model atmospheres and 
the methods to determine the observables will be presented in 
Sect. |3] The approach used to determine the stellar parameters 
(temperature, C/O, mass, log(g)) is presented in Sect. |4] while 
in Sect. |5] we summarise the comparison with the evolution- 
ary tracks. We present a detailed discussion of the single ob- 
jects studied in this work in Sect. |6] followed by a more general 
discussion on the obtained results (Sect. |7]i and the conclusions 
(Sect.EJ. 

2. Observations and data reduction 

Five mildly variable stars were selected for this investigation: 
CR Gem, HK Lyr, RV Mon, Z Psc, DR Ser. In Table[T]the coor- 
dinates of the objects, near-infrared photometry, and the llfim 
IRAS flux are presented together with period (?) and ampli- 
tude (A) of variability. Unless otherwise stated A and P are 



derived from the Gen eral Catalogue of Variable Stars (GCVS; 
ISamus-4 12007 -200^). The amplitude values list the maximum 
amplitude recorded, and might be based only on few data points, 
sometimes recorded on photographic plates. More recent ob- 
serva tions from the public surveys suc h as ASAS (jPoimanskil 
l2002h . and Hipparcos Variability Annex (iBeichman et al.L 19881) . 
show smaller V amplitude variations, often quite stable over 
a long period. These more recent values are listed in the dis- 
cussion on the individual targets (Sect. |6]l. The objects selected 
are semiregular or irregular variables. The variability amplitude, 
and the period indicate that the atmospheres of our targets can 
be represented with hydrostatic models. The literature values of 
A and P are small in comparison to dynamic objects such as 
Mira stars (P longer than one year, and A of several V magni- 
tudes). In the recent grid of dynamic models for C-stars from 
iMattsson et alJ (1201 Ol) almost all the models with a period lower 
than 250 days do not develop a stellar wind. The few exceptions 
are models with extremely low temperature, or very high C/O 
value. The same applies to the prediction s of the previous gen- 
erations of dynamic models (i.e. Table 1 i Gautschv-Loidl et al.L 
2004). The temperature-density structure of these windless mod- 
els is only slightly varying around the hydrostatic configuration. 
These small variations of the atmospheric structure are also re- 
flected in the re sulting spectra,but as can be seen in Fig. A. 1 of 
iNowotnv et al.l (1201 0). the spectra of objects with very mild pul- 
sations and no mass loss do not differ significantly from the cor- 
responding one based on a hydrostatic model. The hydrostatic 
spherical symmetric approximation will be checked with inter- 
ferometric observations (where possible) also in Sect. 12.21 We 
collected for all the targets IJHKL infrared spectroscopy, and K 
broad band interferometric observations. 



2.1. Spectroscopy 

The s pectra for our t a rgets were obtained with the UIST instru- 
ment (iRamsav et all l2004l) on the United Kingdom InfraRed 
Telescope (UKIRT) as part of the UKIRT Service Programme: 
u/serv/1790 and u/serv/1818. The first run of observations 
was on 29 June 2008; the second one on 20, 23, 25, and 26 
January 2009. Four giism spectra were collected for every tar- 
get, coveting the following spectral ranges: IJ with wavelengths 
in the intervaQ [0.862, 1.418] jum; HK with wavelengths in the 
interval [1.395,2.506] fim; shortL, and longL bands covering 
[2.905,3.638] and [3.620,4.232] //m, respectively. We obtained 
also spectra for three standard stars to be used for the calibration 
procedure. The spectral type, photometry, and effective temper- 
ature of the standard stars are listed in Table |2] Since our ob- 
jects were bright, we observed with a 2-pixel (0.24 arcsec) wide 
slit. Targets and standard stars were nodded along the slit, with 
spectra being taken in an A-B-B-A sequence to facilitate first- 
order sky subtraction. Exposure times of 2-coadds x 2 seconds 
or 1-coadd x 4 seconds were used with each giism, depending 
on the brightness of the source. The nodding sequence was ex- 
ecuted twice, resulting in a total integration time of 32 seconds 
per source in each grism. The resolution of the spectra ranges 
between 400 (IJ) and 1800 (longL). HIP 92946 was used to 
remove telluric features for targets observed in 2008 (HK Lyr, 



Although the wavelengths range of the // grism starts at 0.86 fim, 
only the part beyond 0.9 pm is usable because the spectral blocking fil- 
ter blocks emission (essentially has zero transmission) below this wave- 
length and, in any case, the instrumental throughput drops considerably 
towards the / band. 
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Table 1. Targets list, available photometry and variability information. 



ID 


RA 


DEC 


J 


H 


K 


L 


f[12(/m] 


Var. Type 


Amplitude 


Period 








[mag] 


[mag] 


[mag] 


[mag] 


[Jy] 




[mag] 


[days] 


CR Gem 


06:34:23.92 


+ 16:04:30.30 


3.36 


2.07 


1.46 


0.88" 


39.31 


Lb 


1.20 S 


250 


HKLyr 


18:42:50.00 


+36:57:30.89 


3.23 


2.15 


1.62 


0.99* 


22.74 


Lb 


0.40 Vf 


186^ 


RVMon 


06:58:21.49 


+06:10:01.50 


3.06 


1.96 


1.43 


0.85" 


31.13 


SRb 


2.19 B 


131 


ZPsc 


01:16:05.03 


+06:10:01.50 


2.11 


1.10 


0.70 


0.4^ 


33.42 


SRb 


1.30 p 


144 


DR Ser 


18 47 21.02 


+05 27 18.60 


4.12 


2.79 


2.08 


1.38" 


15.9 


Lb 


0.50 V« 


196S 



Notes. The amplitude and the period (unless otherwise stated) are derived from the General Catalogue of Variable Stars (GCVS; Samus 

• photogr a 

■ '"' lEpchtein eTarHl990h : ''" iBagnulol ( 119961) ; '^' iBeichmanet al'Ul988l) ;"'^' lPoimanskil f200j) 



b007-2009 ). The band used for ampb tud e of variability measu re ment is also repo rte d; ' p' stays for photogr aph ic plates . Kerschba um et al, 
jl996gb ; '*> lKerschbaum et"aLHl996bh " ^ 



Table 2. Spectral type, photometry and effective temperature of 
the standard stars used for the spectroscopic calibration. 



ID 


Sp. Type 


V 


J 


H 


K 


Tuff 






[mag] 


[mag] 


[mag] 


[mag] 


[K] 


HIP 92946 


A5V 


4.62 


4.160 


4.163 


4.085 


7 880" 


HIP 33024 


FOV 


5.753 


5.308 


5.232 


5.088 


7 400* 


HIP 5544 


FOV 


5.160 


4.764 


4.545 


4.393 


7 500" 



Notes. T he photometry is from the 2MASS A ll Sky catalogue 
Cutri et al .. 2003); <"' Erspamer & NorthI (l2003l) : <*' [Masana et aL 
20061) ; '"' ISolano & Femlevl '( ll997h . 



DR Ser), while HIP 33024 and HIP 5544 were used for the 2009 
observations (CR Gem, RV Mon and Z Psc, respectively). 

Data reduction was perfonr ied using the ORAC-DR UIST 
pipeline jCavanagh et all l2003h . The subtraction of the sky re- 
moves the OH lines and the thermal background emission. At 
this stage the spectrum is made of a positive and a negative pait 
which can be extracted. An Argon lamp is used for the calibra- 
tion of the object in wavelength. Before every target observation 
a telluric standard was observed. The observations of target and 
reference star are caiTied out always in a range of airmass be- 
tween LO and L3. Only in one case we had airmass L5. The 
airmass difference between target and standard is usually of the 
order of 0.1 thus producing a negligible error in the flux cali- 
bration. Once that the object and the standard were calibrated in 
wavelength, the telluric correction was performed. For this pur- 
pose we used the routine IRFLUX from the StarlinkQ software 
package. IRFLUX divides the standard spectrum by a black body 
with assigned temperature (the one of the standard). The out- 
come of this operation is a spectrum including only the telluric 
features, atmospheric and instrumental effects. Finally the tar- 
get is divided by this last determined spectrum and the telluric 
calibration is performed. The standard stars are hot enough to 
be approximated with a black body in the covered wavelength 
range. 

Once the telluiic connection was performed a few emission 
lines appeared in the spectrum of the targets. Those are resid- 
ual of the atmosphere of the standards which in this part of the 
spectrum does not contribute with continuum (this is flat), but it 
shows few hydrogen lines. We identified those lines by compar- 
ing the spectrum o f the standards with th e corresponding hydro- 
static models from IShulvak et aP (l2004l) . The region where the 
lines are located was not taken into account in the fitting proce- 
dure, as well as the data at the edge of the atmospheric windows 
(shaded region in Fig.[T]i 





1,6 1-8 2.0 2,2 2,4 
Wavelength [fim] 



longL 



5.0 6.2 J.4 J,6 

Wavelength [^m] 



HK Lyr 
CR Gem 
DR Ser 
RV Mon 



J, 7 J, 8 5.2 4,0 4,1 

Wavelength [;.m] 



^ Starlink software is available for downloading on the webpage 
http : // starlink . j ach . hawaii . edu/. 



Fig. 1. Comparison between all the spectra of our targets. The 
main molecular features are indicated. The shadowed (yellow) 
regions are affected by absorption of the Earth atmosphere. 



The reddening correction is always a very delicate task 
in spectroscopy. Every spec trum was dered dened using the 
EXTIN . PRO IDL procedure (lAmores & Lepfnel l2005l) . The code 
computes the visible extinction (Ay) along a path from the 
Sun to any point in the Galaxy, specified by galactic coordi- 
nates and distance. Unfortunately the distance measurements 
associated with our objects are very uncertain. We selected 
three reference distances for each target from t he literature: the 
distan ce measured by the Hipparcos satellite (iPerrvman et al.L 
1 1997b . a distance est i mated assuming a constant K mag- 
nitude dC laussen et al. ', 1987'), and a distance estimated by 
iBer^eat & Chevallier (2005) which is based on Hipparcos mea- 
surements b ut taking into account three distinct biases (see 
Sect. 2.3 of i Bergeat et all l2002h . Hereafter the different dis- 
tances will be indicate with the abbreviation liBer, i^cia, i^Hipp- 

The distances for each object are listed in Table |4] The val- 
ues are in most of the cases very different between each other, 
but no trend or systematic behaviour can be identified. For the 
Hipparcos measurements the eiTors associated to the distances 
are also given. No error estimation was found in literature for 
t/Berg and t/ciau therefore we assume at least 20% of uncertainty. 
This corresponds to 40% error on the luminosity. 

Although the distances have very different values the result- 
ing reddening correction is very similar (as expected since these 
are nearby objects). We tested the effect of the different redden- 
ing correction on the estimation of the parameters (in particular 
the temperature). The typical shift for the temperature determi- 
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one measurement available. The panels in the first row, identi- 
fied with the letter 'a' represent the measurements for DR Ser. 
In the row 'b' the points for HK Lyr are plotted, and in the row 
'c' there are the points for RV Mon. In column ' 1' we present the 
uv coverage. In column '2' the visibility points are presented ver- 
sus the baseline. HK Lyr and RV Mon are observed always with 
the same baseline setup. The variation in the projected baseline 
is due to the Earth rotation. The UD diameter is plotted versus 
position angle in column '3'. In column '4' the UD diameter are 
presented versus the visual phase in order to investigate even- 
tual variations in the angular size which might be due among 
other reason to the pulsation of the stars. The phase estimation 
for RV Mon and DR Ser is based on th e visual light curve from 
the ASAS database dPoimanskil |2002|) : while for HK Lyr it is 
based on the light curve from the Hipparcos Variability Annex 
(Beichman et al., 1988). Small trends can be identified in the 
case of RV Mon and HK Lyr, they may be due to cycle-to-cycle 
variation, but we can also not exclude completely an instrument 
issue (see detailed discussion in Sect.|j 



3. Hydrostatic models and synthetic observables 

For the determination of the parameters we use the new grid 
of hydro static and sp herically symmetric model atmospheres of 
lAringer et al.l (120091) . They are compu ted with COMARCS, a 
progr am based on the MARCS code (iGustafsson et al.L Il975l 



nations is of the order of 10 K, which is low compared to the 
error bar. 

Except for DR Ser, where no correction was applied (see de- 
tailed discussion in Sect. 16.5b . we used for our computation the 
reddening correction estimated for the Hipparcos distance. 

2.2. Interferometry 

Observations for this investigation were primarily taken with the 
Palomar Testbed Interferometer (PTI). PTI was an 85 to 110 
m baselines H and K band (1.6 fim and 2.2 /um) interferom- 
eter located at the Palomar Observatory in San Diego County , 
California, and is described in detail in I Colavita et al] (|1999() . 
It had three 40-cm apertures used in pairwise combination for 
measurements of stellar fringe visibility on sources that range 
in angular size from 0.05 to 5.0 milliarcseconds, being able to 
resolve individual sources > l.O mas in size. PTI had been in 
nightly operation since 1997 and was decommissioned in 2009, 
with minimum downtime during the intervening years. In addi- 
tion to the carbon stars observed as part of this investigation, 
appropriate calibration sources w ere observed as well and can 
be found in Ivan Belle et all ( |2008|) . 

The calibration of the squared visibility (V^) data is 
performed by estimating the interferometer system visibility 
(VgYg) using the calibration sources with model angular di- 
ameters and then normalising the raw carbon star visibility by 
V^Ys estimate t he measured by an ideal interferome- 
ter at that epo ch (Mozurke wich et al.L 119911: iBoden et all Il998t 
Ivan Belle & van Belle . JOOj). Uncertainties in the system vis- 
ibility and the calibrated target visibility are inferred from in- 
ternal scatter among the data i n an o bservation using standard 
error-propagation calculations (iBode n et al. . 1999). Calibrating 
our point-like calibration objects against each other produced 
no evidence of systematics, with all objects delivering reduced 

y2 = i. 

PTI's limiting night-to-night measurement error is ctv^ ~ 

SYS 

1.5 - 1.8%, the source of which is most likely a combination of 
effects: uncharacterized atmospheric seeing (in particular, scin- 
tillation), detector noise, and other instrumental effects. This 
measurement error li mit is an empirically established floor from 
the previous study of iBoden et al.l (1 19991) . 

In our sample there is one star, Z Psc, without PTI obser- 
vations. One visibility measurement taken with IOTA (Infrared 
Optical Telescope Array) is available in literature for this object. 
The visibility point value was obtained in the K filt er (A = 2.2 
um, A A - 0.4 fim) and it is published in Table 1 of iDvck et alJ 
([1996^. We will use this value for our investigation. 

Table lA.ll gives an overview on the observational parameters 
of the interferometric observations used, namely UT date, base- 
line, position angle, calibrated visibility with associated error 

From the relationship between visibility and uniform disk 
(UD) angular size, = [2Mx)/x]\ where x = nB6uDA-\ we iist7.' while' c'^Hr^ndC rdJ0rgens fTsli as OS data. 



2008) in the version used by lj0rgensen et al.l (Il992l) . and 



Aringer et al.l (1 19971) . The models are generated assuming hy 



drostatic local thermal and chemical equilibrium. The molecular 
and atomic opacity is treated in the opacity sampling (OS) ap- 
proximation. The parameters that characterise a model are: ef- 
fective temperature, metallicity, log(§), mass, and C/O. 

3.1. Synthetic spectra 

Our observed spectra have higher resolution than the ones of the 
grid published in Aringer et al. (2009), therefore new synthetic 
spectra were computed for a subset of the initial grid of hydro- 
static models. 

We assumed solar metallicity for our targets. The subset of 
new spectra covers the following parameters 0: 2 400 < T^g < 
4000 K with steps of lOOK; Z/Z© = 1; -1.0 <log(g[cms-2]) < 
-hO.O; M = 1 or2M0; C /O^ 1.05 1.10 1.40. 

The COMA code jAringeii I2OOOI) was used to compute 
the opacities for the different layers of a given temperature- 
density atmospheric structure. The opacity sources included 
in the calculations for th e continuum are listed in Table 1 of 
iLederer & Aringer! (l2009b . Voigt profiles were used for atomic 
lines, and Doppler profiles for the molecules. All the main 
molecular opacities typical fo r C- stars were include d: CO 
("Goorvitch & ChackerianVl994), C2 ('Oue rci et all [19741) . HCN 
(Harris et al.. .2006.) . CN (j0rgensen. 1997t in the form of line 



may establish uniform disk angular sizes for the carbon stars 
observed by PTI since the accompanying parameters (projected 
telescope-to-telescope separation, or baseline, B and wavelength 
of observation A) are well-characterised during the observation. 
This uniform disk angular size will be connected to a more phys- 
ical limb darkened angular size in Sect. |4] By plotting the UD di- 
ameter versus time, we could check the reliability of our approx- 
imation with hydrostatic models. From the UD diameter versus 
position angle the assumption of spherical symmetry is checked. 
In Fig.|2]we present the interferometric observations for three of 
our five targets. The missing objects (CR Gem, Z Psc) have only 



This is in agreement with the previous works (L oidl et al.l 1200 ll 
Aringer et al., 2009; Lederer& Aringer, 2009) with the excep- 
tion of the C2 linelist which was not scaled in the infrared range. 
The scali ng of the s f value s for C2 was suggested for the 



iLoidl 



ime b y IJorgi 
etal.l(l2001h 



first time b y J0rgensen (1 19971) and afterwards introduced by 



for fitting spectra of C-rich stars in the range 
of 0.5 - 2.5 yum. The authors keep the original linelist up to 
1.15 fim, they scale it by a factor of 0.1 beyond 1.5 fim, and use 



^ In the model grid the lower limit of lo g(,^) varies according to the 
temperature. For more details see Table 1 in lAringer et al.l (|2 0091) . 
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DR Ser 



150 




-150 



-150 



150 



50 



-50 



-150 



150 



-50 50 
u[m] 

HK Lyr 



(lb) ^ 



-50 50 
u[m] 

RV Mon 



-150 



-150 



-50 50 
u[m] 



150 




150 




50 100 
Bosellne [m] 



50 100 
Bosellne [m] 



150 




150 




31 48 
PA[degree] 




56 59 
PA[degree] 




Fig. 2. Checking of the hydrostatic and spherical symmetric approximation for the interferometric data. The series of panels ' a' is 
referred to DR Ser data, panels 'b' to HK Lyr, and panels 'c' to RV Mon. The plots in the column '1' represents the Mv-coverage of 
the observations, in column '2' the visibility points acquired are plotted versus the baselines. In column '3' the correspondent UD 
diameter is plotted versus position angle. In column '4' the UD diameter is plotted versus the visual phase. 



a linear extrapolation in the transition region. Dif ferent authors 
poin ted out already the need of new C2 line data. lAringer et al.l 
(l2009i) showed that the C2 scaling does not affect so much the 
model structure, while it introduces a variation in the spectral 
range between 1.3 and 2.1 fim. In the same work a discrepancy 
was observed when comparing synthetic and observed (H - K) 
colours. This discrepancy was investigated in terms of C2 opac- 
ity. The authors show in their Fig. 15 how the photometry ob- 
tained forjmscaled models is in better agreement with observa- 
tions. lAringer et al.l (l2009h conclude that this could be an indi- 
cation that the use of the linelist in its original version is more 
appropriate, as long as there are no other more precise sources 
available. 

In order to check this conclusion further, we performed 
the test as e xplained in Garci'a-Hernandez et al. (2007), and 
lUttenthaler & Lebzelterl (120101) . In most of the cases we ob- 
ta ined a lower result of the test for the Ca-unscaled linelist 
of lOuerci et al ] (fT97l . This result was also judged by eye in or- 
der to test the method. In Fig. [3] the effect of the C2 scaling in 
the HK band is shown. The full line is the spectrum of the star 
Z Psc (upper panel) and HK Lyr (lower panel) convolved to a 
resolution of 200. The dotted line is the synthetic spectrum cor- 



responding to a model with scaled C2 opacity, while the dashed 
line is the same model with the original C2 list. The result of 
the test, t he check by ey e, together with arguments of previ- 
ous work of lAringer et al.l (12009) led us to opt for the original 
C2-unscaled linelist. 

The resulting opacities are used as input for a spherically 
symmetric radiative transfer code which gives the synthetic 
spectra with a resolution of 18 000 in the wavelength range 
0.8 - 25 yum. We convolved these spectra to get same resolution 
as the observed data. 



3.2. Synthetic visibility profiies 

The spherical radiative transfer code produces an additional out- 
put beside the total spectrum: monochromatic spatial intensity 
profile. Synthetic visibility profiles for the /T-broad band PTI 
setup were computed for a subset of the grid of hydrostatic mod- 
els. The metallicity, the temperature range and the C/O ratio 
were fixed (more details in Sect. I4.lt . To compute the visibil- 
ity profiles, i.e. the interferometric observables, we proceeded as 
follows. 
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Fig. 3. Comparison between the observed spectrum of two tar- 
gets (black full line, Z Psc in the upper panel and HK Lyr in the 
lower), a spectrum computed assuming scaled C2 opacity (dotted 
line) and the spectrum with unsealed C2 (dashed line). Shadowed 
bands mark regions of poor atmospheric transmission. 



We defined a set of narrow-band filters centred on the sam- 
pled wavelength of the transfer function for the PTI /T-broad 
band setup. Then the monochromatic intensity profiles were con- 
volved with the filters so defined. The visibility profiles were 
computed as Hankel transform of the narrow -band intensity pro- 
files. The spatial frequencies were converted from AU"' to base- 
lines in meters, for a given distance. The squared visibility broad 
band profile finally is obtained as: 

= - CI 'I 



were the sum was calculated over the narrow band filters; S , 
is the transfer function of PTI in the K broad band setup; F/ is 
the flux integrated over the narrow band filter, and finally V, is 
the visibility corresponding to the narrow band filters. Following 
this approach for the computation of the broad band visibil- 
ity profil e, the bandwidth smearing effect is properly taken in 
accou nt dKervella et all 120031: IVerhoelstl \200^ iPaladini et all 
I2009h . 



4. Parameter determination 

All the spectra of the targets appear very similar (Fig. [1]. They 
are dominated by the CN, C2 and CO bands between 0.9 and 
2.5 yum, while the main absorption feature at 3 fj.m in the L-band 
is due to HCN and C2H2. In Fig.[T]the regions affected by telluric 
absorption are shadowed. Thes e regions were of c ourse not con- 
sidered in our fitting procedure. iLoidl et al.l(l2001l) computed the 
flux ratio of various hydrostatic models (e.g. Fig. 3 of their pa- 
per) to investigate the effect of changes in the stellar parameters. 
They demonstrated that s pectral feature s are r nainly sensible to 
Teif and C/O. Moreover, lAringer et aP (l2009l) showed that the 
synthetic spectrum below 2.5 fim of stars with T^g > 3 000 K 
mildly changes for different values of \og{g) and mass. 

Therefore we determine Tgff and C/O from spectroscopy, and 
we turn to the PTI interferometric observations to determine 
the remaining characteristic parameters of the models: mass and 
log(^). 

4. 1. Temperature and C/O ratio 

The ratio between the 3 jum (C2H2) and the 5.1 ;um (C3) fea- 
tures is the most powerful tool for d etermining temperatur e 
and C/O ratio for hydrostatic C-stars dj0rgensen et al.l l2000l) . 
Unfortunately our spectral coverage does not reach the C3 fea- 
ture at 5.1 /im. 

Following the approach of L oidl et al.l (12001) we use the en- 
ergy distribution of the available spectra as main temperature 
indicator The observed spectrum was convolved to a lower res- 
olution of R = 200. Every target of our list was compared with 
a grid of synthetic spectra with the same resolution. The grid 
of models is dense enough in temperature to allow an accurate 
statistic approach which is described in detail in Appendix lAl 
The Tgff value was estimated for every wavelength range sep- 
arately, except for the longL band. This range of the spectrum 
is mostly flat and free from features therefore not very sensitive 
for the temperature estimation. We computed also a temperature 
taking into account all the four parts of the spectrum. The re- 
sulting temperatures are presented in Table |3] Figure shows 
the mean squared deviatiorQ derived by fitting first single pieces 
of the spectrum ('77", "HK", "shortL") and then by fitting all 
the pieces of the spectrum at once ("all the wavelength range") 
for HK Lyr Our intention is to look for possible minima in the 
temperature distribution. 

The temperature values obtained from the different ranges of 
wavelength show a trend: by using the IJH bands as indicators, 
cooler temperatures are obtained. This trend can be hardly ex- 
plained with problems during the calibration. The data reduction 
procedure (see Sect. 12.1b was accurate enough that we do not 
expect, in the overall flux distribution, systematic errors which 
are large enough to explain the systematic variations in the tem- 
perature. The temperature difference might be explained as the 
effect of an optically thin dust shell surrounding the star. The ef- 
fect of a dusty shell on the ene rgy distribution of a star is shown 
in Fig. 3,4 of iNowotnv et al.l (1201 ll) . The model used for this 
simulation correspond to a very dynamic star. Our objects are 
rather static, but we speculate that a thin dust shell would affect 
the spectrum in the same way with scaled-down intensity. This 

"* One root mean square value (rms) was derived by fitting each 
model. The models were divided in bins of temperature; the single ms 
values were squared (mean square deviation) in order to be summed for 
obtaining one value for every bin. The number of models in every bin is 
not constant, therefore the final mean square deviation was normalised 
by the number of models in every bin. 
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Fig. 4. The four plots show the mean squared deviation obtained 
by fitting single portions of the spectra (IJ, HK, shortL), and 
by fitting all the portions of spectra at once (all the wavelength 
range) for HK Lyr. These plots demonstrate the accuracy of the 
shortL measurement. 



shell would absorb the light in the IJH band and emit it at longer 
wavelength. As a direct consequence a temperature determina- 
tion based on IJH band is systematically lower as can be seen in 
the upper panels of Fig.|4] If the shell is thin enough the effect on 
the short edge of the L band is negligible. This makes the shoitL 
band temperature determination quite robust, as can be seen in 
the lower left panel of Fig. |4] 

We note that this is not the case if the star is a strongly 
pulsating variable (such as Mira). In this case the profile of the 
3 i^m feature, which is the result of the su perposition of molecu- 
lar op acity and dust continuous emission dGautschv-Loidl et al.L 
|2004|) . it will be sensitive to the dynamic processes of the 
atmosphere. The respective temperature derived for every star 
by fitting the shortL band will be assumed hereafter 

The second quantity which mainly modifies the appearance 
of a C-rich spectrum is the C/O ratio. Unfortunately the rela- 
tively low resolution of the spectra available limits the precision 
of this measurement. The CO band at 2.29 yum, the C2 bands at 
1.02 and 1.20/im, and the CaHaH-HCN at 3 n m were con s idered 
as indicators for the C/O ratio. Following to iLoidl et al.l (1200 Ih 
the C2 features longward of 1 .20 fim were not used because of 
the uncertainty in the C2 opacity. Two reference wavelengths 
were chosen to indicate start and end of every selected band. In 
low resolut ion spectroscopy no continuum window is available 
for C-stars (Palad ini et al.Ll2009l) . therefore a pseudo-continuum 
was obtained by linear interpolation of the points at the edge 
of the chosen band. The selected portion of spectrum was nor- 
malised to this pseudocontinuum. Following this approach we 
derived the equivalent width for the observed spectrum and for 
all models with the temperature determined from the previous 
step. Grouping the models in sets according to the C/O (i.e. every 
"C/O group" includes models with same temperature, but differ- 
ent \og(g) and mass), a root mean square value was obtained for 
every C/O by comparing the observed equivalent width and the 
synthetic one. 

In the cols. 6, 7, 8 of Table [3] the resulting rms for the three 
C/O values (1.05, 1.10, and 1.40) is given. The minimum rms is 
evidenced in bold. All the stars of the sample have a low value of 



C/O ratio which corresponds to 1.05, except for HK Lyr where 
the C/O = 1.4. 



4.2. Mass, log(g) and distance 

As already pointed out by other authors it is hard to detect the 
effect of mass and \og{g) from low resolution spectroscopy of 
C-rich stars. This is confirmed in our seiies of plots in Fig. l6lfT0l 
where models with the same parameters except mass and \og(g) 
are compared with the spectra of the targets. The temperature 
and C/O ratio of the models were determined following the pro- 
cedure previously described. 

Due to these considerations we decided to treat the observ- 
ables of interferometry as completely independent quantities. As 
interferometry aims to measure the radius of the target, the dis- 
tance becomes an important parameter Unfortunately the dis- 
tance measurements available for these objects are rather impre- 
cise (typical error of the order of 20%) and often in contradiction 
with each other with differences between measurement relative 
to the same objects which are larger than eiTors. The problem we 
face is degenerate since we have to deal with radius, distance, 
mass, and all these quantities are related to each other 

We handle the problem in the following way. Three distances 
obtained with different methods were chosen from the literature 
for every object (see Sect. 12.11 and Table |4), and a set of syn- 
thetic visibility profiles was scaled to every distance. The set has 
a fixed temperature, and C/O ratio. We obtained for every dis- 
tance the combination of stellar parameters (M, log(g)) from the 
best fitting profile. 



Z = 0.019 




3.60 



3.55 



3.50 



3.45 



3.40 



3.35 



3.30 



Fig. 5. Zoom into the H-R Diagram where C-stars are located. 
The solar metallicity isochrones (Maiigo et al., 2008) are plot- 
ted in grey and the small numbers indicate the predicted present 
mass. The sampling of points in the isochrones is shown with 
tiny squares. The position of the C-stars is indicated by a shad- 
owed area, the uncertainty about this area is discussed in Sect.|5] 
The different symbols plotted correspond to the parameters de- 
termined for every star and distance assumed: plus for CR Gem, 
diamond for HK Lyr, star for RV Mon, triangle for Z Psc, and 
square for DR Ser 
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Table 3. Summary of temperature, C/O determination, and rms values coiTesponding to the different C/O. 



ID 






T,s(HK) 


reff(shortL) 


r^ff (all spec) 


rms 1.05 


rmsi.io 


nnsi.4o 




[K] 




[K] 


[K] 


[K] 








CR Gem 


2 700 ± 


200 


2 860 + 200 


3 070 ±50 


2 920 ± 190 


0.021 


0.027 


0.038 


HKLyr 


2740 + 


220 


2 920 ±210 


3 090 ±50 


3 080 ± 120 


0.015 


0.016 


0.013 


RVMon 


2950 ± 


300 


2 930 ± 230 


3 170 ± 50 


3 210 ± 140 


0.006 


0.008 


0.009 


ZPsc 


3 000 ± 


330 


3 080 ±210 


3 130 ±60 


3 170 ± 130 


0.014 


0.017 


0.016 


DR no redd correct 


2 790 ± 


250 


2 820 ± 200 


3 080 ±40 


3 030 ± 170 


0.011 


0.013 


0.014 



Notes. The bold face corresponds to the best values of temperature and C/O. More details in Sect. 14. II 

Table 4. Distances measurements and stellar parameters of the model which best fits spectroscopic and interferometric measure- 
ments. 



Distance 


CR Gem 


HKLyr 


RVMon 


ZPsc 


DR Ser 




920 pc 


730 pc 


670 pc 


465 pc 


1295 pc 




L = 1 1 000 L© 


L = 7 186 L© 


L = 6400 L© 


L = 4534L© 


L = 18 050 L© 




Teff = 3 100 K 


Teff = 3 100 K 


T,g = 3 200 K 


Teff = 3 100 K 


Teff = 3 100 K 




C/O = 1.05 


C/O = 1.4 


C/O = 1.05 


C/0= 1.05 


C/0= 1.05 




M = 2Mq 


M = 2Mq 


M = IMq 


M = 2Mq 


M = 2Mq 




log(g) = -0.40 


log(g) = -0.20 


logis) = -0.40 


log(g) = +0.0 


log(g) = -0.60 




780 pc 


900 pc 


1000 pc 


600 pc 


990 pc 




L = 9025 Lq 


L= 11 389Lo 


L= 12 932Lo 


L = 7 170 Lq 


L = 1 1 390 L© 




Teff = 3 100 K 


Teff = 3 100 K 


T,g = 3 200K 


Teff = 3 100 K 


Teff = 3 100 K 




C/O = 1.05 


C/O = 1.40 


C/0= 1.05 


C/0= 1.05 


C/0= 1.05 




M= IMq 


M = 2Mq 


M = 2Mo 


M = IMq 


M = 2Mq 




log(g) = -0.60 


\og{g) = -0.40 


\og(g) = -0.40 


log(g) = -0.50 


log(g) = -0.40 




323!f;; pc 


1 369!^'; pc 


450!;:;; pc 


323!-^ pc 


1 690+-'*'" 








L = 2 574 Lq 


L = 2267 Lq 










r,ff = 3 200 K 


T,tf = 3 100 K 










C/0= 1.05 


C/0= 1.05 










M = IM© 


M = IMq 










logfe) = +0.0 


log(g) = +0.0 





Notes. The bold face distances correspond to the most likely values found by our stellar parameter analysis. 



5. Comparison with evolutionary tracl<s 

At this point of our investigation we have one combination of 
stellar parameters for every object and every distance. To con- 
strain the choice of the parameters further, we compare them 
with the recen t isochrones for thermally pulsing AGB stars from 
iMarigo et al.l (2008). The sele ction of the isochrones follows the 
same criteria as lMarigo et al.l (1200 8): ages log(t/yr) between 7.8 
and 10.2, and the spacing in log(f) is 0.1 dex. 

Fig. |5] shows a zoom into the region of the H-R diagram 
where AGB stars are located. The shadowed area identify the 
region where C-stars are expected. Due to the discrete sampling 
of the points in the isochrones (small square), the transition re- 
gion from M- to C-stars is not very well defined. According 
to the stellar evolution calculations (P. Marigo, private commu- 
nications) all the stars on the cool side of the "hook" of the 
isochrones are C-stars. Some values of the present mass are 
marked on the isochrones: they will be compared with the mass 
of the best fitting hydrostatic models (see Sect.|7]i. 

We overplot to the isochrones the points corresponding to 
the parameters we determined. Every star is indicated with a dif- 
ferent symbol centred on the temperature and luminosity associ- 
ated to the model: a plus for CR Gem, an asterisk for HK Lyr, 
diamond for RV Mon, triangle for Z Psc, and square for DR Ser 
The error bar of the temperature are the ones associated to the 
shortL band Tgff determination, therefore they are centred on the 
Tetf (shortL) values given in col. 3 of Table |3] The errors on the 
luminosity are not plotted here to avoid confusion but they are 



of the order of 40% (corresponding to the uncertainty on the dis- 
tance measurements liBei-g, i^ciau)- 

6. Discussion on individual targets 

6.1. CR Gem 

CR Gem is classified as an irregular Lb variable with a variabil- 
ity amplitude in the B band of 1.20 niagint he General Catalogue 
of Variable star s (GCV S; ISamus-tl 12007 -2009). Nevertheless 
IWhitelock et al.l (l2008h classify this star as SRb with a period 
of 250'' a nd the ASAS V band light curve shows an amplitude of 
0.5 mag (PoimanskiL 120021) . The spectral classification given in 
the GCVS isC8,3e(N). 

The 3 jum shortL temperature we deter mine is around lOOK 
higher than the 2 960K measured by iBergeat & Chevallieil 
(2005). . The C/O we get is 1.05, we could not find other es- 
timates for this value in literature. The observed spectrum is pre- 
sented in Fig. |6] It is in general well reproduced with a few ex- 
ceptions: the range 1.2 - 1.4/im and the /T-band between 2.1 
and 2.3 jim where the real data show an excess compared to the 
models (see general discussion in Sect.]?). 

For this star we collected only one point of PTI visibil- 
ity (Fig. |6] lower right panel), therefore it is not possible to 
check for asymmetries. All the visibility profiles coiTespond- 
ing to the Hipparcos distance are not extended enough to fit 
the visibility point, therefore this distance can be excluded. For 
the Claussen distance the best fitting model has mass 1 Mq and 
\og{g) = -0.60, the luminosity is 9 000Lo ™d '^^e radius of 
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the star 330 /?0. For the Bergeat distance, the largest value of 
distance estimated for this star, the best fitting model has 2 Mq 
and \og(g) - -0.4, the luminosity is 1 1 000 Lq and the radius 
310 Rq. 

6.2. HKLyr 

According to the GCVS the visual amplitude of HK Lyr is 
1.80 mag, the spectral classification C6,4(N4), and the star 
is cl assified as irregular Lb. The Hipparcos Variability Annex 
jBeichm an et al., 1988) reports for this star a period of 196'' 
and an amplitude of variability of 0.4 in the Hipparcos magni- 
tude. From optical spectroscopic analysis this star is enr i ched i n 
Lithium and Technetium (Abia et al., 2002; Boffin et ~anil993h . 

From our spectroscopic investigation the temperature and 
C/O ratio obtained are higher than the ones estimated 
in literature (T a ble |3] l. While our terr iperature is 3 080 K, 
lOhnaka & Ts^ (1199^ give 2 866 K and iBergeat & ChevalHeii 
(2005) es timate 294 5 K. Concerning the C/O ratio we obtain 
1.4 while lAbia et al.1 (12002) estimate 1.02. The fit of the spec- 
troscopic data is quite successful except in the region between 
1.2 and 1.45 //m (Fig. |7] upper left panel). 

Eight visibility points are available from PTI observations. 
The UD diameter computed by fitting the single points do not 
differ from each other. Over the time interval and position angle 
of the observations, the star does not show notable variation nor 
evidence of asymmetries (Fig.|2] row 'b'). We fitted all the visi- 
bility points with synthetic profiles scaled for different distances 
(Fig Q lower right panel) . 

None of the models corresponding to the Hipparcos distance 
can fit the points, therefore we excluded this distance. The best 
fitting parameters corresponding to (ifierg are 2 Mq and \og(g) - 
-0.2. This model has a luminosity of L = 7 200 Lq and a radius 
R = 300 Rq. The best fitting model corresponding to c/ciau has 
again mass 2 Mq but \og(g) = -0.4. This model has a luminosity 
of L = 1 1 400Lo and a radius R = 370 Rq. In Fig.|5]the resulting 
parameters for HK Lyr are plotted with asterisks. 



the UD versus position angle. A careful check was performed 
to exclude any problem coming from the calibration procedure, 
but still we cannot exclude completely that the fluctuation in the 
visibility is due to an instrument issue. From the astrophysical 
point of view the observed trend could be explained in differ- 
ent ways: (i) it could be an asymmetry plus a temporal variation 
(for example a large convective shell), (ii) or an effect of the pul- 
sation (panel '4c'). The temporal variation due to the primary 
short period of variability is not enough to explain the jump in 
the visibility. We derived from the light curve a difference in 
phase for the primary period of 0.16 between the two sets of ob- 
servations. Accor ding to the predictions of the dynamic model 
atmospheres (Pala dini et al.L i2009: see Fig. 7 and Fig. 9,) a hot 
C-star should show a difference in UD-diam of 0.06 mas, defi- 
nitely smaller than 0. 16. The two sets of observations have been 
obtained with a time difference of 3.5 times the length of the 
long secondary period. This long term variation may, thus, be 
responsible for the observed variation in size. Further observa- 
tions are necessary to reach a conclusive interpretation. Keeping 
this in mind we performed a fit of the interferometric points with 
a hydrostatic model atmosphere as a first step. This might be 
followed by further investigations. From the lower right panel 
of Fig. [8] it is immediately obvious that no synthetic visibility 
profile can fit at the same time all the observations. The best 
fitting model obtained for the Hipparcos or Claussen distances, 
can fit only the single point observed in 1999. The other models 
are either too extended or not extended enough to fit the other 
observations. The best fitting model for dHipp has M = IMq, 
log(g) = +0.0, L = 2 600 Lc^ and ^ 170 Rq. For the distance 
determined by IClaussen et al. l (fT987l) the single data point is fit- 
ted by the model wifli 2Mq, log(g) = -hO.4, L = 13 000 Lg and 
a radius R - 370 j?f7i. On the o t her ha nd for the distance deter- 
mined bv lBergeat & Chevallierl (l2005b the three points observed 
in 2008 are fitted by the model with M - IMq, log(g) - -0.4, 
L - 6400Lo, and R - 260 Rq. Nevertheless the effect of vari- 
ability or possible asymmetries need to be investigated in more 
detail, and the parameters derived for this object should be con- 
sidered with caution. 



6.3. RVMon 

RV Mon is a C4,4-C6,2(NB/R9) star, its variability class is SRb 
with a primary pulsation period of 131'' and a long secondary 
period detected by Houk ( 1963) of 1047"'. The amplitude of vari- 
ability given in the GCVS is 2.19 mag in the B-band, while the 
ASAS light curve suggests a V amplitude of 0.3 mag. 

The temperature obtained with our fitting procedure is in 
between the values given in the literature. lOhnaka & Tsuiil 
(IT996I) get Teff ^ 3 33 0K, we obtained 3 170K while 
iBergeat & Chevalliej (l2005h give 2 9 1 K. No estimation of C/O 
was found in the literature to be compared with our 1.05. The 
model we adopted for RV Mon matches quite well the spectro- 
scopic observations in all the ranges except in the region of the 
CO bands (longward 2.29 yum) where the observations show an 
IR excess compared to the synthetic spectra. The intensity of the 
C2 feature at 1 .77 yum is not well reproduced but this is also at 
the border of an atmospheric window. 

For this star we collected 4 visibility points. The first one 
is in 1999, the other three in 2008. Although the observations 
were carried out at very similar baseline, the visibility jumps by 
0. 1 between the two epochs of observation. This jump is visible 
also in the UD-diameter where it amounts to ~ 0.16 mas which 
is larger than the estimated error on the UD diameters ( ~ 0.06 
in the worst case). This difference is also noted in the plot of 



6.4. Z Psc 

Z Psc is an SRb variable with amplitude derived from photo- 
graphic plates of 1.3 mag. The period in the GCVS is 144''. The 
ASAS light curve has an amplitude of 0.52mag i n V band. T he 
spectral c lassification from the GCVS (ISamus4l 12007-20091) is 
C7.2(N0). lAbia et al.l(l2002l) . and lBoffin et al.l(ll993h detemiiried 
a '^C/'^C= 55, moreover they found enhancement of Tc, and no 
trace of Li in the spectrum. This implies that Z Psc is a standard 
low-mass TP-AGB star, and an intrinsic C-rich star 

M any temperature estim ates were given f or this star in liter- 
ature : lLamber t et al.'('l986) obtained 2870 K.'Dvck et al.' ('I99r' 



3240 K, lOhnaka & Tsuii ( 1996) 3150 K, Bergeat & Chevallie^ 
(l200l 3095 K. All these values were obtained with different 
methods. Our measurement (see Table |3) is very close to the one 
given bv lOhnaka & Tsuiil (Il996h The C/O = 1 05 we estimated 
is in agreement with the 1 .014 of'Lambert et al. '(fT986^, and with 
the more recent 1.01 value of Abia et al. (2010). The match be- 
tween models and observations is shown in Fig. |9] From spec- 
troscopy a small disagreement is barely detectable in the C2 band 
at 1 .20 /im, while the 3 yum feature of the model is too deep to fit 
the data (the star is hotter than the model). 

Z Psc is the only one among our targets with no PTI data. 
For t his star we used a n observation available from the litera- 
ture dDvck et all 1 19961) taken in the /T-broad band with IOTA 
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Fig. 6. Comparison of the UKIRT/UIST spectra (black line) and PTl interferometric measurements (black triangle) of CR Gem 
with hydrostatic models predictions. The upper left panel illustrate the IJ range of the spectrum, the upper right panel shows the 
HK wavelength range. The lower left panel shows the shortL range. The shadowed bands mark the region with poor atmospheric 
transmission. The molecules that contribute to the spectrum are also indicated. The synthetic spectra are overplotted in grey (blue, 
and red in colour version). The lower right panel shows the interferometric data point with the models which best fit the data for 
different distances overplotted. 



(2.2 yum) as already described in Sect. 12.21 The fit of the in- 
terferometric IOTA data is presented in the lower right panel 
of Fig. |9] The best fitting model for the Hipparcos distance 
has a mass M = IMq, log(§) = -hO.O, L = 2 300LfD, and 
R - 1 65j^(T). The distance determined bv iBergeat & ChevaUierl 
(l2Q05h gives a best fitting model with the same log(g) obtained 
for Hipparcos distance fit but M - 2Mq, L - 4500Lo ^^'^ 
a radius R - 230 Rq. In the case of the Claussen distance 
we have the best fitting which is the model with mass 1 Mq, 
log(g) = -0.5, L = 7 150 Lq, and 7? = 300 Z?©. 



6.5. DR Ser 

According to the GCVS, DR Ser is an irregular variable Lb with 
an amplitude in the B-band of 2.99 mag. The ASAS light curve 
of this object is quite stable over the last 2000 days, it has a 
period of 196'', and a V band amplitude of 0.3 mag. The spectral 
class is C6,4(N). 



This object is the one in our sample with the highest value 
of interstellar reddening. Applying the reddening correction we 
noted that the observed spectrum looks completely unrealistic. 
Remembering that reddening correction for single objects is al- 
ways very uncertain we decided to skip the step of the reddening 
correction for DR Ser Further investigation and a detailed map 
of the ISM in the region of this star are needed. 

We determined a temperature of 3080 K. In the lit- 
erature the Teff values associated to this object ar e usu - 
ally lo wer (see Sect. Hi: 2 500 K (lAbia & IsernL I 997h 



2570 K (lThompsonet alJ.l2002h. 2650 K l lSchoier etall l2005l: 
iBerpeat & Chevameiil2005l) . lAbia & Isernl (ll997l) estimated the 
isotopic ratio '^C/'^C = 6, which is also surprisingly low! The 
C/O we fou nd, 1.05 is lower than the C/O - 1.26 obtained in 
hterature bv lEglitis & Eglitel (Il995b . 

There are 54 PTI points available for this object in a period 
between June 2000 and August 2001. The UD diameter deter- 
mined agree with each other within a range of +0.02 mas, there- 
fore we considered for the model fitting all the points together 
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(Fig. [Tol l. None of the profiles computed at the Hipparcos dis- 
tance could fit the data points. The best fitting models have a 
mass M - 2 Mq, and log(g) = -0.4 or -0.6 for dckm and dserg 
respectively. The first model has L = 11 400 L0 and a radius of 
370 «o; the second one L = 1 8 000 L© and 7? =470/?©. It is im- 
pressive to note how all the visibility points observed at different 
epochs and position angles fits the hydrostatic synthetic profile. 



7. Discussion 

Based on the fitting procedure previously described, we can 
identify a sequence in effective temperature for our targets. The 
temperature increases starting from CR Gem, DR Ser, HK Lyr, 
Z Psc, RV Mon. Our temperature determination is very precise 
thanks to the large grid of models available and to the use of 
the shortL band spectrum. The Teff here determined are always 
compatible within the error given in literature. 

The C/O ratio determination is not as accurate as the one of 
the temperature because of the very coarse spacing in the model 
grid, but mainly because of the low resolution of the spectra. 
Although the C/O values we determine have to be considered as 
indicative, they are quite in agreement with the ones estimated 
in literature. We obtained the same C/0= 1 .05 for all the stars 
except HK Lyr which has C/0= 1 .4. Those values are in agree- 



ment with the range [1.01 , 1.41 reported in hter ature for classi- 
cal C-stars (see Fig. 42 of 'L ambert et all (Il986l) . and Table 2 of 
Bergeat & ChevalHer (2005)). 

In general there is a good agreement between synthetic and 
observed spectra. Small discrepancies are observed mainly in the 
region between 1.2 and 1.4 fim, between 2.1 and 2.3 fim and on 
the left edge of the 3 fim feature. They can be attributed to small 
calibration problems, but also to a small amount of dust. This 
last point is supported by the fact that the differences detected 
in the IJ band between 1.2 and 1.4 fim are larger for the lower 
temperature stars (CR Gem, HK Lyr, DR Ser). We simulated 
the effect on the s pectrum of a sma ll amount of dust by using 
the DUSTY code dlvezic & Eliztuit Il997l) . Indeed the IJ band 
of the spectrum fits slightly better for stars which have a large 
J - K value (CR Gem, DR Ser). We also checked the effect of 
this small amount of dust on the stellar parameter determination. 
The effect on the temperature and on the interferometric observ- 
able was within the error bar we already took into account, there- 
fore we proceed with the hydrostatic models. A detailed analysis 
of spectro-interferometric observations with hydrostatic+dusty 
model atmosphere will be presented in a follow-up work. 

From the fit of the interferometric observations with mod- 
els at different distances we note that in three over five cases 
no model can fit the observations at the distance estimated from 
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Hipparcos. This is not such a surprise considering the error as- 
sociated to these distance measurements. Further consideration 
about the distance can be made by studying the position of the 
best fitting models in the H-R diagram (Fig|5]). 

Some of our hydrostatic models lie in a region of the H-R 
where higher masses are expected. If we consider the agreement 
between the two masses (from hydrostatic models and evolu- 
tionary calculations) as a criterion to select the correct distance 
of the objects then we have the following: 

- CR Gem. The point with lower luminosity is the one fit- 
ting the range of mass, therefore the distance determined by 
Claussen can be retained, and the one of Bergeat discarded. 

- HK Lyr Both the models have the same mass and lie in a 
region in between 1 and 2.6 Mq. Considering an error on 
both the distances of 20% on the distances (i.e. a factor 1.44 
in luminosity) none of the values can be excluded. 

- RV Mon. The Hipparcos distance is in the correct region 
of mass, but too far from the C-stars region, therefore can 
be discarded. Both the other two distances correspond to a 
model which fits the region of masses, therefore they can be 
both retained. 

- Z Psc. In this case the Bergeat distance can be discarded 
since the best fitting model has 2 Mq , but falls in the 
range of 1 Mq. Considering the error on the luminosity both 



the distances from Hipparcos and Claussen can be retained. 
Nevertheless the position of the model at the Claussen dis- 
tance is more probable since this is closer to the C-stars re- 
gion. 

- DR Ser The distance estimated from Bergeat can be dis- 
carded since the model hes in the region of the stars with 
3Mq. 

The parameters with the most probable distances, with asso- 
ciated errors in temperature and luminosity are shown in Fig . [TTI 

Taking into account the uncertainty associated to the tran- 
sition phase M-C stars (see Sect. |5]l, all the stars have at least 
one combination of parameters that lie in the region where the 
C-stars sequence starts. This is in agreement with the hydro- 
static scenario for our targets: the objects are supposed to be hot 
C-stars, with relatively low mass-loss and low C/O. They "re- 
cently" turned into C-stars and while moving on the evolution- 
ary track they will increase the production of carbon and will 
lose more and more mass because of the stellar wind. Our work 
is partially limited by the fact that only model atmospheres with 
mass 1 and 2 are available. 

Two objects of our sample show some peculiarity which de- 
serve further investigations. The interferometric observations of 
RV Mon show a possible non-hydrostatic nature. All the visi- 
bility points acquired with PTl cannot be fitted with the same 
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model. This might be due to some instrumental effect, while the 
physical explanation would be an evidence of an asymmetry or 
(more probably) of a variability effect. According to the interfer- 
ometric measurements, the star DR Ser is an hydrostatic object. 
We were able to fit with the same model 54 visibility points ob- 
served in different epochs and position angles. On the other hand 
the fit of the IJ and HK band is not very good, and the reddening 
correction of the spectrum makes the fit even worse. A detailed 
ISM map for this region of the sky is needed. 



8. Conclusions 

More than a simple parameters determination, this work aims 
to identify which of the stellar parameters of C-stars can be 
constrained by using different techniques like spectroscopy and 
interferometry. Which limits appear by investigating only with 
one tool? Infrared interferometry is a relatively young field of 
research, but it can already give a significant contribution on the 
stellar parameters determinations. 

In this work we confirm that low resolution infrared spec- 
troscopy is not sensitive to \og(g) and mass determinations of 
C-stars. We show clearly that a highly precise determination of 
the temperature can be achieved for C-stars with no significant 



dust contribution thanks to the use of the 3 fim feature (C2H2 -1- 
HCN). The determinations based only on the IJH bands tend to 
underestimate the temperature. 

The hydrostatic model atmospheres adopted in this work 
were able for the first time to fit simultaneously spectro- 
scopic and interferometric observations. This is a very impor- 
tant achievement considering that inaccuracies in model atmo- 
sphere reflect negatively in the check of the stellar evolution 
calculations. On the other hand the inaccuracy in the distance 
determination for this class of objects reflects negatively on the 
parameter determination, as already noted in Wittko wski et al.l 
(l2Q01h . We considered a criterion to select the distance on the 
agreement between mass of the hydrostatic models and mass 
predicted by the isochrones. Ac cording to this, the distance esti- 
mated bv lClaussen et al.l (Il987h is the most probable among the 
three distances we used for this work. Taking into account also 
the uncertainty in the transition phase from O-rich to C-rich star 
in the evolution models, we obtained that all the objects of our 
sample have at least one combination of observationally deter- 
mined parameters which fits the C-stars region. We were able to 
associate temperature, C/O, mass, log(g) and a range of reason- 
able distances for all the objects of our sample. 

A grid of models with larger coverage of the mass range (at 
least from 1 to 4 Mq), as well as a more dense grid of points 
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for the isochrones, would improve significantly our estimations. 
The advent of a new dedicated mission (Gaia) might improve the 
distance problem for C-stars. 

The synergy between different techniques of investigation 
and realistic theoretical atmospheric and stellar evolutionary 
models will be the key to determine precise characteristics for 
this class of objects, and to be able to understand clearly the 
physical processes which drives their evolution. 
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Appendix A: Statistical approach for temperature 
determination 

In order to assign an observed spectrum to one out of a sample of 
available models, a statistical approach is employed: each model 
represents a Gaussian-shaped sandheap the width of which (in 
terms of standard deviation) is determined by the point-to-point 
scatter p of the spectrum. The height of the sandheaps can be ad- 
justed by statistical weights. Thus the shape of the nth sandheap 
is given by the Gaussian probability density function 



4 



(A.l) 



The "distance" d„ between the observed spectrum and the model 
is the rms deviation between the observed and the synthetic spec- 
trum. The statistical analogy of the model assignment is to pick 
out a single grain of sand at the given distances from the syn- 
thetic spectra. Now the question is, "What is the probability of 
the picked grain to belong to the «th sandheap?" 

The total number of available grains at the selected position 
is determined by the sum of densities for all models. 



y.Gidn) 



(A.2) 



which represents our normalization condition. Thus the relative 
amount of grains belonging to model n - and hence the probabil- 
ity of th e observed spectrum to represent this model - evaluates 
to 

Gid„) 



Pr„ 



(A.3) 



In the context of our application, the only parameter of in- 
terest is the effective temperature, whence the temperature in- 
terval covered by the models is divided into a set of intervals 
I,, t - 1, ...,T. Since the number of models belonging to such a 
temperature bin is not unique, the marginal probability density 
of temperature has to be corrected for this bias. In the sandheap 
metaphor, we would have to provide each temperature interval 
to be represented by a unique number of grains, and the proba- 
bility of our spectrum to be assigned to the temperature interval 
/, evaluates to the weighted sum of individual model assignment 
probabilities 



Pr(/,) = 



(A.4) 



Er=iPr*(//) 

m, and denoting the number of models belonging to the interval 

This normalisation is demonstrably valid, because the inte- 
gral probability for all temperatures evaluates to 



(A.5) 



List of Objects 

'CR Gem' on pagep 
'HK Lyr' on pageQ 
'RV Mon' on pagelj 
'Z Psc' on page|2 
'DR Ser' on page|2 
'HIP 92946' on page 
'HIP 33024' on page 
'HIP 5544' on page[3| 
'Z Psc' on page 
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Table A.l. Interferometric observations. 



ID 


UT Date 


B 

[m] 


PA 
[deg] 


V ±crv 


Additional informations 


CR Gem 


11/03/1999 


101.49 


80 


0.324 ± 0.004 




HKLyr 


05/20/1999 


109.79 


71 


0.277 + 0.004 






05/21/1999 


109.78 


78 


A TT7 ■ A AAC 

0.277 ± 0.005 






05/22/2000 


109.29 


63 


0.258 ± 0.005 






05/22/2000 


109.33 


64 


0.265 + 0.004 






05/22/2000 


109.76 


71 


0.238 ± 0.010 






(Jj/zZ/ZUUU 


1 An nn 
109. / / 


/2 


A O /I /I 1 A AA'J 

U.Z44 + U.UU3 






iJj/ZZ/ZiJUiJ 


1 no Tj 




A o/io _i_ n nn'2 






()3/22/z()()() 


109.72 


8 1 


0.242 ± O.OOj 




RV Mon 


1 1 /A /I / 1 AAA 

11/04/1999 


100.37 


63 


f\ A A C 1 AA1C 

0.446 + 0.015 






A A C A A O 

09/25/2008 


108.68 


52 


A OAC A AAC 

0.305 + 0.005 






A A C 1^ A A O 

09/25/2008 


106.98 


55 


A OOA 1 A A A /I 

0.320 + 0.004 






A A 1'^ C /'^ AAO 

09/25/2008 


104.79 


57 


A O /I A ■ A f\r\£i 

0.349 ± 0.006 




ZPsc 


10/08/1995 


38.21 




0.810 ± 0.123 


(1) 


DRSer 


06/03/2000 


106.29 


55 


0.387 + 0.008 






06/19/2000 


107.05 


54 


0.535 + 0.024 


(2) 




06/19/2000 


106.01 


56 


0.497 ± 0.035 


(2) 




06/19/2000 


105.45 


56 


0.487 ± 0.032 


(2) 




06/19/2000 


104.83 


57 


0.515 + 0.024 


(2) 




06/20/2000 


103.33 


59 


0.435 + 0.013 






06/20/2000 


102.19 


60 


0.429 ±0.013 






04/18/2001 


104.41 


57 


0.407 ± 0.010 






04/18/2001 


104.20 


58 


0.415 + 0.004 






04/18/2001 


103.23 


59 


0.407 + 0.003 






04/18/2001 


103.06 


59 


0.413 + 0.007 






04/18/2001 


102.13 


60 


0.410 ± 0.005 






04/23/2001 


106.14 


55 


0.418 ± 0.013 






04/23/2001 


103.56 


58 


0.410 + 0.012 






04/23/2001 


100.63 


62 


0.429 + 0.011 






04/23/2001 


99.18 


64 


0.435 + 0.011 






04/23/2001 


97.73 


66 


0.440 ± 0.012 






04/29/2001 


105.27 


56 


0.389 + 0.006 






04/29/2001 


105.01 


57 


0.391 + 0.005 






04/29/2001 


103.67 


58 


0.396 + 0.005 






04/29/2001 


102.84 


59 


0.406 + 0.005 






04/29/2001 


102.68 


59 


0.407 ± 0.008 






05/03/2001 


86.43 


15 


0.587 + 0.014 






05/03/2001 


86.40 


15 


0.597 + 0.015 






05/03/2001 


85.92 


15 


0.591 + 0.019 






05/03/2001 


84.84 


16 


0.593 + 0.018 






05/03/2001 


83.03 


17 


0.626 ± 0.023 






05/04/2001 


86.07 


15 


0.583 + 0.009 






05/04/2001 


85.98 


15 


0.580 + 0.009 






05/04/2001 


84.76 


16 


0.599 ± 0.010 






05/04/2001 


82.50 


17 


0.615 ± 0.007 






05/04/2001 


82.24 


17 


0.621 + 0.011 






05/05/2001 


86.31 


15 


0.589 + 0.006 






05/05/2001 


85.79 


15 


0.598 + 0.007 






05/05/2001 


84.85 


16 


0.605 ± 0.005 






05/05/2001 


83.54 


17 


0.613 ± 0.008 






05/05/2001 


81.86 


17 


0.626 + 0.006 






05/23/2001 


104.92 


57 


0.402 + 0.007 






05/23/2001 


1 AO OA 

103.80 


58 


A /I 1 ^ 1 A A A /I 

0.416 + 0.004 






06/21/2001 


103.71 


58 


0.404 ± 0.006 






06/21/2001 


102.06 


60 


0.420 ± 0.005 






06/21/2001 


100.30 


63 


0.433 + 0.005 






06/21/2001 


98.60 


65 


0.460 + 0.009 






07/16/2001 


79.72 


14 


0.627 + 0.012 






07/22/2001 


105.69 


56 


0.402 ± 0.019 





Paladini et al.: Determination of stellar parameters of C-rich hydrostatic stars 

Table A.l. continuated. 



ID UT Date 


B 

[m] 




PA 
[deg] 


V + crv 




Additional informations 


07/22/2001 


105 


.55 


56 


0.405 + 


0, 


,019 


07/22/2001 


104, 


.73 


57 


0.393 + 


0, 


,015 


07/22/2001 


104, 


.59 


57 


0.413 + 


0, 


,016 


07/22/2001 


101, 


.55 


61 


0.427 + 


0, 


,016 


07/22/2001 


100, 


,54 


62 


0.438 + 


0, 


,014 


08/22/2001 


86, 


,22 


15 


0.586 + 


0, 


,022 


08/22/2001 


86, 


,04 


15 


0.581 + 


0, 


,020 


08/22/2001 


85, 


,95 


15 


0.596 + 


0, 


,022 


08/22/2001 


85, 


.41 


16 


0.586 + 


0, 


,026 


08/99/9001 






1 f> 




V. 


090 


08/22/2001 


80 


.33 


18 


0.631 + 


0, 


,028 


08/22/2001 


79, 


.97 


18 


0.661 + 


0, 


,055 


08/23/2001 


98 


.03 


66 


0.456 + 


0, 


,011 


08/27/2001 


86 


.28 


14 


0.564 + 


0, 


,005 


08/27/2001 


86 


.45 


14 


0.565 + 


0, 


,003 


08/27/2001 


86 


.46 


15 


0.566 + 


0, 


,004 



a^ lDvcketai]n99a) 

(2) spurious values, discarded of the computation 
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